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ABSTRACT 

We present a method of estimating perturbative coefficients in Quantum Field 
Theory using Pade Approximants. We test this method on various known QCD 
results, and find that the method works very well. 
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By using the first n coefficients in a series expansion, we have estimated the 
(n + l)-st perturbative coefficient in Quantum Field Theory (QFT). Though there 
is currently no theoretical basis for extrapolating coefficients in the perturbative 
loop expansion of QFT by our method, our results have thus far been in good 
agreement with the calculated coefficients of quantum electrodynamics (QED), as 
well as with series in statistical physics, condensed matter theory and mathematics 
[1-6]. In this paper we compare our method to the perturbative loop expansion of 
quantum chromodynamics (QCD) at the five-loop level. We shall present results 
for the i?-ratio, the R T ratio, the QCD j3 Function and two QCD Sum Rules. 

Our method makes use of Pade Approximants (PA) and enables us to obtain 
an estimate and an error-bar for each coefficient. We define the PA 



a Q + a l X + -.- + a N X N 

[ 1 1 i + blX + ... + b M xM [L) 



to the series S where we set 

oo 

[N/M] = S + 0(X N+M+1 ) and S = S nX n . (2) 

One solves Eq. (2) and then predicts the coefficient of the next term Sjv+M+i ■ 
This is what we do in this paper. One can also use the full PA to estimate the sum 
of the whole series S. This is what we will do in the future [7,8]. For a detailed 
description see Refs. [5] and [6]. 

The PA's are known to accelerate the convergence of many series by including 
the effects of higher (unknown) terms, thus providing a more accurate estimate 
of the series. We have recently proved the following theorem, which provides a 



useful sufficient condition for the PA's to be accurate. Defining f(n) = £nS n and 
considering 

9(n) = (3) 
a sufficient condition for the PA's to converge is that 



lim g(n) —0 . (4) 

The PA's thus provide reliable estimates of asymptotic series whose coefficients 
diverge as 

S n = n \ , (5) 

as is believed to be the case in QED and QCD [9]. It can easily be shown that 
Eq. (4) is satisfied for S n given by Eq. (5). In the cases of these and other series 
whose Borel transform has a finite radius of convergence, the higher-order PA's 
give progressively better approximations to the Principal Value of the transform 
integral over Borel singularities. It is easy to check that for series with one or two 
simple Borel poles, (i.e. IR and UV renormalons) , Pade approximants predict the 
next term in a given series with a rapidly increasing precision. In these cases an 
analytic estimate ~ M\/N M can be made for the relative error of the [N/M] Pade 
approximant prediction of the next term in the series. 

It might be objected, however, that these cases are not sufficiently complicated 
to be realistic. Therefore, as an exercise, we have evaluated Pade approximants 
to the large- Nf limit of the vacuum polarization D-function in QCD, which is 
known to all orders in a s , and whose Borel transform contains an infinite series of 
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double poles at both positive and negative integers. 10 Once again, as seen in the 
Figure, the Pade approximants' predictions of the next term in the series converge 
rapidly, in agreement with the above-mentioned estimated error for the [N/M] 
Pade approximant. The convergence of the Borel transform of the D-function 
series, in particular, indicates that our PA approach is well suited for perturbation 
series with the asymptotic behavior expected in QCD. 
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Relative errors in the [N/M] Pade approximants: (a) to the QCD vacuum polarization D-function, 
evaluated to all orders in the large- Nf approximation 10 - the rate of convergence agrees with 
expectations for a series with a discrete set of Borel poles, and (b) to the Borel transform of the 
D-function series, where the convergence is particularly striking. 



4 



We now turn to some QCD applications of our PA approach. Let us first 
consider the R T ratio 11 where R T is defined as follows: 

R T = r(r — > v + hadrons)/r(r — > ez/z/) = 3SEw( r r + r r) (6) 

where Sew = 1-019 is the electroweak correction and r\ = —1.58% is the non- 
perturbative contribution. The perturbative QCD contribution is given by 

r> = 1 + — + 5.202 f — ) 2 + 26.36 (-Y + (109.2 ± 12.9) ( (7) 

71 \ 7T / V 7T / \ 7T / 

where the last coefficient is our estimate of the five-loop contribution. Using the 

12 13 

experimental average ' value R T = 3.623(17), we obtain for the strong coupling 
constant 

a s (M T ) = 0.325(6) . (8) 

where the error does not include various systematic uncertainties, which go be- 

g 

yond the scope of this letter and are discussed elsewhere. Numerically Eq. (7) 
becomes r T = 1 + 0.1035 + 0.0557 + 0.0292 + 0.0125=1.201(19) . One can see that 
the perturbative series is converging, albeit somewhat slowly. Using the known 
(3 function (see Table II for our estimate of the four-loop (3 function), Eq. (8) 
corresponds to A^ 3 ) = 355(11) MeV in the MS Scheme. Stepping up through the 
4 and 5 fermion-thresholds (m c = 1.5 GeV and m& = 5 GeV, respectively) we get 
A( 4 ) = 306(11) MeV and A( 5 ) = 218(9) MeV, and, hence, 



a s (34 GeV) = 0.1399(11) 
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(9) 



For the R ratio, we have 
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R = 3SQ 2 f r : r = 1+ (—) +1.409 



V 7T / 




-12.805 



( ^_ (8 , 5±10 ,) 

(10) 



The last term in Eq. (10) is our estimate for the five-loop contribution to R where 
we here extrapolated the related Adler D-function. Using Eq. (9) we obtain 

r(34 GeV) = 1 + 0.0445 + 0.0028 - 0.0011 - 0.0003 = 1.0459(4) . (11) 

This series for r in Eq. (11), where the contributions in each order, up to five loops, 
are given, seems to converge nicely. Experimentally there are two measurements of 
r(34 GeV). They are r(34 GeV) = 1.049(7) 15 and r(34 GeV) = 1.056(8) 16 . It can 
be seen that the extrapolation prediction in Eq. (11) is in good agreement with 
these experimental values. We can now evolve a s up to M z , the Z boson mass. 
Our result is 



which is consistent with the latest experimental value obtained from total cross- 

17 

section measurements at LEP 



For r at M z we get r(M z ) = 1 + 0.0378 + 0.0020 - 0.0007 - 0.0002 = 1.0389(2) . 

We now present our estimates for higher-order perturbative coefficients for R, 
R T , the QCD f3 Function and two QCD Sum Rules. All of the results presented 
here are in the MS scheme. The corresponding results for other schemes can also 



a s (M z ) = 0.119(2) . 



(12) 



a s (M z ) = 0.126(6) . 



(13) 
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be obtained. However the scheme- dependence should decrease in higher order, 
disappearing if all orders are known. 

To provide systematic-error estimates, we first apply our method to S n , ob- 
taining our estimate We then apply it to the reciprocals r n = 1/S n , and 
take the reciprocal again to obtain S^h. We then take differences t n = r n+ i — r n 

(2) 

and apply our method to obtain a third estimate S^^. Our systematic-error es- 



timate is A/2, where A = 



0(2) _ o(l) 
D n+1 °n+l 



. We then combine the diagonal and 



non-diagonal estimates of S^_ v weighted by 1/A 2 , to obtain our final estimate for 
S n +i- 

In Table I we present our results for R and R T . The first entry in each case 
is the four-loop result. The light-by-light contribution is small but should be 
added to i?(MS). Our estimates of the four- loop coefficients, based on PA's to 
lower-order coefficients, agree well with the known exact results, providing a sound 
phenomenological footing for our method. The second entries are our estimates 
for the five-loop coefficients and the numbers in brackets are our estimate of the 
systematic errors. The results —96.8 KS and 105.5 KS are the estimates of Ref. [18], 
obtained using a completely independent method, Optimized Perturbation Theory 
(OPT). The agreement with our estimates is very good. 

In Table II we present our results for the QCD /3-function. The agreement with 
the known 3-loop results is very good, and we present the first estimates of the 
four-loop QCD (3 function. Note that the three- and four-loop results are scheme- 



dependent, and we use the MS scheme throughout. The 3-loop result is the same 



for any MS-type scheme. Since we use the MS result for the 3-loop coefficient our 



estimate for the four-loop coefficient is also for the MS scheme. 



In Table III we present our results for the Bjorken sum rule for deep inelastic 
(unpolarized) neutrino-nucleon scattering. To compare with Ref. [18], we multiply 
their results BjnSR by —2/3. Our results for O(a^) are very good again, and again 
the agreement with Ref. [18] in O(a^) is excellent. We present our results for the 
Bjorken sum rule for deep inelastic polarized electron-nucleon scattering, — BjpSR 
in Table IV. Again our results in 0(o£.) are good and our 0(af) estimates agree 
with Ref. [18]. The Gross-Llewellyn Smith sum rule (GLSSR) differs from the 
BjpSR by the light-by-light contribution: GLSSR = BjpSR - 0.413 / . The light- 
by-light contribution here and for R should be treated separately. This contribution 
is small, however, for all cases of interest. 

In Table V we present our results for the R ratio for various Nf. The small 
difference for Nf = 5 compared with Table I is due to a slight difference in averaging 
the various estimates. The results for the four-loop coefficients are excellent and 
the agreement with Ref. [18] for the five-loop coefficients is satisfying. 

Although it is interesting that our estimates for the next term agree with 
Ref. [18] and we believe both of us are right, we cannot be certain. Both of us may 
be wrong. Fortunately, even a crude estimate of the next term is sufficient since 
their contributions to the full series are small. However, it is important to know 
they are small! 

In conclusion, we have used our estimation method, which makes use of Pade 
Approximants, to estimate various perturbative coefficients in QCD. Our estimates 
for the known terms is very good. Moreover our estimates for the next unknown 
terms agree very well with the results of Ref. [18] in all cases where comparison is 
possible. 
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Note added: 

A phenomenological extraction of the five-loop D-function coefficient from mea- 
sured moments of r decay data has recently appeared 19 , which confirms our pre- 
diction. 

In the process of our analysis we are comparing our results with explicit expres- 
sions given in ref. 14 and in papers listed in ref. 18. The original references for the 
NNLO calculations are: R(s) - Ref. 14 and S.G.Gorishny, A.L.Kataev, S.A.Larin 
in "Standard Model and Beyond: from LEP to UNK and LHC" Dubna, USSR, 
October 1-5,1990, S. Dubnicka, et al., Eds. World Scientific (1991) p.288; the QCD 
beta-function - O.V. Tarasov, A. A. Vladimirov and A.Yu. Zharkov, Phys. Lett. 
B93(1980)429; the Bjorken non-polarized sum rule - S.A. Larin, F.V. Tkachov and 
J. Vermaseren, Phys. Rev. Lett. 66(1991)862 and Bjorken polarized sum rule 
- S.A.Larin and J. Vermaseren, Phys. Lett. B259(1991)345. We thank Andrei 
Kataev for correspondence on this matter. 
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TABLE I 



Estimates for R(MS), N f = 5, and R T {MS), N f = 3, at the four- (first row) 
and five-loop (second row) order. The numbers in brackets are the estimated la 
error-bars. The four-loop results are compared with the exact (known) results and 
the five-loop results are compared with those of KS. Nf is the number of fermions 
(quarks). 



SERIES 


ESTIMATE 


EXACT 


R(MS), N f = 5 


-10.20(1.53) 


-12.76 




-87.5(10.8) 


-96.8 (KS) 


R T (MS), N f = 3 


27.06(6.77) 


26.37 




109.2(12.9) 


105.5 (KS) 



TABLE II 

Estimates for the QCD /^-function in three- (first row) and four-loop (second row) 
order. The 3-loop results are compared with the exact (known) results. Nf is the 
number of fermions (quarks). 

SERIES ESTIMATE EXACT 



QCD (3 FCN 

N f = 3 -455(228) -644 

-5920(1956) 



N f = 4 -316(158) -406 

-3058(875) 



N f = 5 -195(49) -181 

-845(105) 
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TABLE III 

Estimates for the Bjorken sum rule for deep inelastic (unpolarized) neutrino- 
nucleon scattering. The 0{a s [) (first row) results are compared with the exact 
(known) results and the 0{a s v ) (second row) results are compared with those of 
KS. We multiply their results by —2/3 for ease of comparison. 



SERIES ESTIMATE EXACT 
BjnSR X - 2/3 

Nf = 3 -13.0(6.0) -18.6 

-116(69) -133 (KS) 

N f = 4 -10.5(4.7) -13.4 

-67.2(17.7) -75.8 (KS) 

Nf = 5 -8.3(3.6) -8.5 

-30.6(1.4) -29.4 (KS) 



TABLE IV 

Estimates for the Bjorken sum rule for deep inelastic polarized electron-nucleon 
scattering. The 0{a s v ) (first row) results are compared with the exact (known) 
results and the 0{a s v ) (second row) results are compared with those of KS. Their 
result should be multiplied by —1 to compare with ours. 

SERIES ESTIMATE EXACT 



- BjpSR 





= 3 


-12.8(6.4) 


-20.2 






-112(33) 


-130 (KS) 


N f 


= 4 


-10.6(4.4) 


-13.9 






-58.2(15.2) 


-68.1 (KS) 




= 5 


-8.5(3.4) 


-7.8 






-21.1(3.4) 


-17.8 (KS) 
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TABLE V 



Estimates for i?(MS) for various Nf at the four- (first row) and five-loop (second 
row) order. The four-loop results are compared with the exact (known) results and 
the five-loop results are compared with those of KS. 

SERIES ESTIMATE EXACT 



R RATIO 



N f 


= 3 


-14.1 ±2.0 


-10.27 




-119.3 ±6.3 


-128.4 (KS) 


N f 


= 4 


-12.2 ± 1.4 


-11.5 




-115.2 ±5.4 


-111.8 (KS) 


N f 


= 5 


-10.07 ± 1.99 


-12.76 




-86.5 ± 10.8 


-96.8 (KS) 
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